Magnetic Flux Periodic Response of Nano-perforated Ultrathin Superconducting Films 
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We have patterned a hexagonal array of nano-scale holes into a series of ultrathin, superconducting Bi/Sb films 
with transition temperatures 2.65 K < Tco <5 K. These regular perforations give the films a phase-sensitive 
periodic response to an applied magnetic field. By measuring this response in their resistive transitions, R{T), 
we are able to distinguish regimes in which fluctuations of the amplitude, both the amplitude and phase, and 
the phase of the superconducting order parameter dominate the transport. The portion of R{T) dominated by 
amplitude fluctuations is larger in lower Tco films and thus, grows with proximity to the superconductor to 
insulator transition. 



The superconductors originally considered by BCS exhib- 
ited spectacularly sharp transitions from a finite resistance to 
zero resistance as a function of temperature 1 1]. Fluctuation 
effects were negligible. Presently, a great deal of attention 
focuses on low superfluid density superconductors for which 
fluctuations strongly influence and substantially broaden their 
phase transitions. These include the high temperature super- 
conductors 01, and in particular, their underdoped versions 
and ultrathin superconducting films near the super- 
conductor to insulator transition (SIT) |6, 7, .8., 9, IQl. For the 
latter, resistive transitions, R{T), can develop widths com- 
parable to or greater than the apparent mean field transition 
temperature, Tco LI 1.1 ■ 

To discuss the effects of fluctuations on the R{T) it is help- 
ful to consider the two component superconductor order pa- 
rameter i/j — |-0o|e"^. In bulk elemental superconductors, the 
sharp R{T) reflect the near simultaneous appearance of a fi- 
nite amplitude, \ipo\, and long range coherence of the phase, 
(f>. In low superfluid density superconductors, however, the 
amplitude first forms at high temperatures and phase coher- 
ence develops at lower temperatures [3 fl4il . High 
on a transition quasiparticles transiently form Cooper pairs, 
which enhances the quasiparticle or fermionic contribution to 
the conductivity, cr/. These pair or amplitude fluctuations give 
rise to the initial drop in R{T) \ 15]. At very low values of R, 
a substantial Cooper pair density exists and the conductivity 
of these bosons, df,, controls R{T). ai, is limited by the mo- 
tion of vortices which causes fluctuations in the phase of the 
cooper pair condensate. In between, a two fluid model may 
best describe the transport fTl. 

Physical interpretations of the R{T) in high sheet resis- 
tance films especially those near the SIT relies heavily on 
distinguishing the amplitude and phase fluctuation dominated 
regimes 1^00 01- Most often, explicit models for (Tf{T) 
and <Jh{T) do not exist as a guide and qualitative arguments 
must prevail. In this paper we present the magnetic flux re- 
sponse of the R (T) of very low superfluid density (high sheet 
resistance) fill llfill films patterned with a nanoscale array 
of holes. We use the quality of the flux response at differ- 
ent points along their transitions to determine the presence of 
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Figure 1: Superconducting transitions for a perforated film. Inset 
shows an SEM image of one of the substrates used in these experi- 
ments. The hole lattice constant is a = lOOnm and Rhoie = 33nm. 



well-defined vortices 1 171 and thus provide insight into where 
the amplitude fluctuation dominated transport gives over to 
phase fluctuation dominated transport on an R{T). In addi- 
tion, we propose that this method can be used to directly de- 
tect the existence of vortices in settings where their presence 
has been contentious 181. 

As a consequence of phase coherence, magnetic field pe- 
riodic behavior is expected in films with multiply connected 
periodic geometries (see inset of Fig. Q 1 18, 191. For exam- 
ple, the resistive transitions of thick Nb films 1 201 on nanopore 
arrays and thick Al films with much larger lattice and hole di- 
mensions 1211 oscillate in temperature while maintaining their 
shape in an applied magnetic field. For holes separated by seg- 
ments on the order of or less 1201 
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(1) 



where ATc is the shift in the transition, is the (dirty) co- 
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herence length, a is the hole to hole spacing ||22|], n is an in- 
teger, and (p is the applied magnetic flux. 0o = Hm ^"^a^^ 
is the magnetic flux corresponding to having one flux quan- 
tum per unit cell and defines the unit cell matching field, Hm- 
Eq. ^is derived within the context of mean field theory, or in 
the absence of appreciable fluctuations of the superconduct- 
ing order parameter and without considering vortex motion. 
Therefore, the physical picture which leads to Eq. ^is that of 
immobile vortices which penetrate the film through the array 
of holes. These vortices generate screening currents which, at 
H = uHm, nearly cancel in the interior of the film so that 
ATc ~ 0. At incommensurate fields the screening currents 
give rise to pair breaking and thus, a reduction in Tc. The 
microscopic picture described by Eq. ^is the familiar Little- 
Parks physics. 

Eq. n accounts well for the flux periodic behavior of thick 
Nb films 1201 . Our data, however, show that this physical 
picture is not valid for ultrathin high sheet resistance films 
with low superfluid density. Amplitude fluctuations prevent 
the formation of screening currents as described above when 
the fraction of the normal state resistance r = R/Rn 0.5. 
In addition, as r decreases (r < 0.15) and screening currents 
appear, the films' susceptibility to thermal fluctuations allows 
activated vortex motion to dominate the transition. Our data 
also indicate that the region of the transition dominated by 
amplitude fluctuations grows with decreasing Tco and hence 
increasing proximity to the SIT. 

Our experiments were conducted on homogeneous quench 
condensed Bi/Sb films similar to those employed in previ- 
ous studies of the SIT and Thermally Activated Flux Flow 
(TAFF) [llj. These films have very low superfluid density 
because they are thin and disordered fill \wk . Two films 
were made simultaneously for each experiment. One film was 
deposited on a fire-polished glass substrate while the other 
was deposited on a nano-perforated Anodic Aluminum Ox- 
ide (AAO) substrate (see the inset of Fig. [Q. The latter as- 
sumed the honeycomb geometry of the substrate. The prepa- 
ration of the AAO substrates can be found elsewhere in the 
literature 26]. Those used in these experiments 

had a hole lattice constant, a, of 100 nm and a hole radius, 
Rhoie — 33 nm. Both substrates were precoated with 50 nm 
of thermally evaporated Ge (in an attempt to smooth the small 
surface roughness of the AAO substrate) before Au contact 
pads were deposited at room temperature. Subsequently both 
substrates were mounted in our cryostat where homogeneous 
Bi films are fabricated at T = 8K by first evaporating a thin 
film of Sb (< Inm, which ensures the films' homogeneous 
morphology) and then depositing the desired thicknesses of 
Bi through sequential depositions. In this way, a series of Bi 
films were fabricated without breaking vacuum or warming. 

The R{T) were measured using standard four-terminal, 
low-frequency ac techniques and acquired in the regime of 
applied currents where the films exhibit an ohmic response. 
Nonlinearities associated with the phase transition to the zero 
resistance state, which are expected at lower temperatures and 
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Figure 2: Iso-dissipative measurements of temperature vs magnetic 
field (normalized by the matching field) for a representative perfo- 
rated film with dsi = l.QOnm, Tco ~ 3.70K, rhoie = 33nm, 
and hole spacing a — lOOnm. Inset shows the corresponding super- 
conducting transition with the arrows indicating the value of R/Rn 
where the curves in the main panel were taken. Eq. Q predicts 
that each curve should exhibit the same magnitude oscillation of 
^ O.OIK. 



have been observed in some wire arrays do not appear in this 
experiment's range l27i E5E9ll . The normal state resistances, 
Rn, of neighboring film regions were compared to assess 
film homogeneity and found to agree to < 7%. Magnetic 
fields were applied perpendicular to the planes of the films. 
Our thermometry consists of a calibrated carbon glass resis- 
tor which has a negligible magnetoresistance in the range of 
fields used in these experiments. 

Data from a series of four superconducting films with Tco 
ranging from 2.65 to 5 K and normal state resistances, Rn = 
R{8K), ranging from 3.3kJl to 300^1 are shown in Fig. 
The reference film exhibits the same range of Tco and a sim- 
ilar increase in the widths of its transitions with decreasing 
Tco- The latter characteristic has been ascribed to a growing 
fluctuation dominated regime [llj. The systematic reduction 
in Tco is believed to result from disorder enhanced Coulomb 
repulsion effects that grow with increasing Rn and possibly 
drive the SIT 130, 31]. 

Despite the breadth of the superconducting transitions they 
do oscillate with increasing magnetic field. The amplitude of 
the oscillations, however, depends strongly on the reduced re- 
sistances, R(T)/Rn, at which they are measured. That is, 
unlike thicker films the shape of the transitions change in field 
1I2O,. .211. Figure|2]shows a typical series of isodissipative (i.e. 
fixed r) measurements of AT |32]. {Rm does not change in 
this range of applied fields.) The lower r curves exhibit os- 
cillations that diminish in size with increasing field and are 
superimposed on a nearly quadratic background. The peaks 
in the data appear at H = nHm — n(2145 ± 50) Gauss 
where n is an integer Hm corresponds closely to the calcu- 
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Figure 3: Normalized amplitude of the first temperature oscillation as 
a function of fractional resistance for 4 perforated films with2.65 < 
Tco < 5.06K, Rhoie ~ 33nm, and a — lOOnm. Inset shows the 
unexpected behavior of the A{r) with T^o near re- 



lated value for one flux quantum per unit cell of the hole array 
(see f2^). Thus, these oscillations result from the collective 
response of the hole array rather than from the responses of 
single holes issll . A feature of particular interest is that the 
peak amplitudes are larger at lower r. 

We characterize the flux periodic response using a normal- 
ized oscillation amplitude, A, for the first oscillation: 



A = 



T{r,0)+T{r,HM) T{r,HM/2) 



2Tco 



Tco 



(2) 



Each solid curve of A as a function of r shown in Fig. |3l 
was calculated using individual measurements of R{T) at 0, 
Hm/2, and Hm- The dashed curve represents the mean field 
oscillation amplitude, Aa/j? — {^o/2a) , derived from Eq. [fl 
5o — 10 nm represents the average value for these films as 
determined from measurements of the upper critical fields of 
the reference films. It varied by 20% over this range of Tcq. 

The oscillation amplitudes exhibit a rich dependence on r 
and Tco, unlike the constant mean field prediction. For all Tcq, 
A grows nearly logarithmically from zero with decreasing r 
with a dependence that is stronger for the films with lower 
TcQ. For r ^ 1, the oscillations surpass 10 times the mean 
field result. The r at which oscillations first appear, ro, is well 
defined and is lower for lower Tcq. Interestingly, the A{r) for 
different Tco cross in the vicinity of a single point near r = 
Tc ~ 0.15. The existence of this crossing point is supported 
by data on another hole size. Its significance is brought out 
partially by the inset in Fig. |3] which shows how A depends on 
Tco at fixed values of r. Vc delineates the portions of the R{T) 
in which A increases or decreases with the primary energy 
scale characterizing the transition, Tcq. 

Accounting for the very large A at low r and the vanishing 
of A for r ~ i clearly requires consideration of the effects 
of fluctuations |34]. Unfortunately, detailed theoretical pre- 



dictions for the R{T) do not exist. Our identification of two 
scales in the data, tq and Tc, however, provide a useful guide 
for the discussion. As we emphasize below, the two scales de- 
fine points on the transitions at which the dissipative processes 
and character of the fluctuation effects qualitatively change. 

Insight into the large amplitude of the oscillations at r < 
Tc comes from investigating the tails of the superconducting 
transitions in magnetic field. Figure |3 shows Arrhenius plots 
of the R{T) collected at H = Q, Hm /2, and Hm for two 
films with different Tco. The tails of the R{T) follow 



R{T) = Roe 



i-To/T) 



(3) 



which is the signature of Thermally Activated Flux Flow 
(TAFF). The energy barrier, Tq, is lower by more than a factor 
of 2 at = i?A//2 than at H = Hm- This nonmonotonic 
variation in To gives rise to the AT oscillations that contin- 
uously grow as r ^ 0. At H = Hm the vortices come 
into registry with the holes and thus their lattice is made stiff 
against thermally activated flux flow. At half integer values of 
H/ Hm, vortices cannot come into registry with the holes and 



10« 



10-' 



10-" 



10-' 



10^ 



1CH 







H=o\ 

T^=2.eSK 




r 1 ! 1 1— 








,\ 

1 1 1 


V 

S 



1.1 



1.2 



1.3 



1.4 



Figure 4: Arrhenius plots for the perforated film with Rhoie = 
33nm. Tco = 2.65/C (top) and Tco = 3.7K (bottom). The slope 
of the curve characterizes the energy barrier against vortex motion 
in the film. The barrier in the perforated film is non-monotonic with 
field. Dashed lines are fits to Eq|3| 
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simultaneously form a stable lattice. Consequently, the acti- 
vation energy is substantially lower at incommensurate fields. 
That lower Tcq films have larger oscillation amplitudes in this 
regime results from the difference in the activation barriers 
for different Tco films as seen in Fig. |3 Even on a normalized 
temperatuere scale, the lower Tcq film has broader transitions 
than the higher Tco film. Correspondingly, the activation en- 
ergy increases faster than linearly with Tcq. This behavior is 
consistent with the condensation energy dictating the activa- 
tion energy scale. A more detailed discussion of the TAFF 
behavior will appear in future work. 

Moving up the transition to r > Vc (Fig. |3}, the oscilla- 
tion amplitudes wane and eventually go to zero. The behav- 
ior as r ^ To can be attributed to enhanced vortex mobility 
that causes them to lose registry with the holes. Nearer Tco 
the shrinkage of the order parameter makes the barrier to in- 
terhole vortex motion smaller and the likelihood of amplitude 
fluctuations greater. Given that the constrictions in the film are 
quasi- Id (a few across) both of these effects will enhance 
interhole vortex motion. These effects are more pronounced 
in lower Tco films since the condensation energy, which is 
proportional to Tco, governs their strength. Consequently, A 
is lower in lower Tco films. As r exceeds rg, all evidence of 
well defined vortices that could be pinned by holes and give 
rise to a periodic response disappears. This last regime can be 
characterized by amplitude fluctuations that become so strong 
that they destroy fluxoid quantization altogether 

Thus, even in the absence of a detailed model of R{T) 
which accounts for fluctuations in these films, we can iden- 
tify three distinct fluctuation regimes defined by rc and tq. 
For r < rc, thermally activated flux motion causes phase fluc- 
tuations. In the intermediate regime, ro > r > rc, the vortex 
mobility is augmented by strong amplitude fluctuations. And 
finally, for < r < 1 amplitude fluctuations dominate and 
vortices do not exist. 

The decrease of ro with Tco implies that the size of the am- 
plitude fluctuation dominated regime increases closer to the 
SIT 1 1 1 1. This trend can be attributed to the growth of both 
the quantum critical regime of the SIT |35| and the classical 
critical regime |36| as Rn increases. Rn ~ i?c/2 for the low- 
est Tco film in the series. It suggests that the fermionic degrees 
of freedom strongly influence the approach to the SIT in uni- 
form films. Experiments on lower Tco, higher Rn films will 
provide insight into whether bosonic (vortices) or fermionic 
degrees of freedom dominate the SIT. 

We propose that techniques similar to those applied here 
may be useful for detecting the existence of vortices in other 
unconventional situations. For example, evidence from Nernst 
Effect 1 3] and magnetoresistance |4] measurements suggests 
the presence of vortices well above Tc in under doped high 
Tc compounds. Underdoped films patterned with an ordered 
array of pinning centers (e.g. holes or magnetic impurities) 
should exhibit flux periodic behavior according to our results. 
Similar arrangements |8| could test the proposal that there are 
vortices on the insulating side of the superconductor to insu- 
lator transition. 



In summary, we have quench condensed nano-perforated, 
homogeneously disordered Bi films. The regular perforations 
induce a phase-sensitive periodic response to an applied mag- 
netic flux and thus provides a probe capable of distinguish- 
ing between the two types of superconducting fluctuations. 
In these low superfluid density films, we have found that the 
strength of this response crosses over between regions of the 
R{T) dominated by fluctuations in the amplitude, of the 
superconducting order parameter to that dominated by fluc- 
tuations in the phase, (p. In addition, as Tco is lowered and 
the system approaches the quantum critical point of the SIT, 
the region of the R{T) dominated by amplitude fluctuations 
grows. 
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